Moving discrete breathers in a Klein— Gordon chain with an impurity 
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We analyze the influence of an impurity in the movement of discrete breathers in Klein-Gordon 
chains. We observe that the moving breather can cross the impurity, can be reflected by it, or can 
be trapped originating a quasi-periodic breather. We flnd that resonance with a nonlinear localised 
mode centred in the impurity is a necessary condition in order to observe the trapping phenomenon, 
as a difference with the resonance condition with a linear localised mode when this problem is 
studied within the Nonlinear Schrodinger Equation approximation. 

PACS numbers: 63.20.Pw 63.20.Ry 63.50.+X 66.90.+r 
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I. INTRODUCTION 

A nonlinear phenomena which has been paid a lot of at- 
tention during recent years concerns to the localisation of 
oscillations in discrete nonlinear Klein-Gordon lattices. 
Since the discovery of intrinsic localised modes or dis- 
crete breathers by Sievers and Takeno |^ , they have been 
found in a great number of systems. MacKay and Aubry 
proved the existence of discrete breathers in Klein- 
Gordon lattices and established a method for calculating 
I them as exact solutions of the dynamical equations (for 
'a review, see ||^, |^). In addition, these localized oscil- 
lations, under certain conditions, can move and trans- 
port energy and they are usually called moving breathers 

JiSii- 

A method used for calculating discrete breathers in 
Klein-Gordon lattices consists in approximating the dy- 
namical equations by a Nonlinear Schrodinger Equation 
, (NLS) and the moving breather by an envelope soliton. 
' One problem that had been studied using the NLS limit 
concerned to the interaction of moving breathers with an 
■ impurity supposing that the breather has a small ampli- 
' tnde However, moving breathers involve oscillations 
,of large amplitude, which implies that the NLS approxi- 
mation may not be as accurate as expected. The aim of 
, this paper is to study the influence of an impurity in the 
movement of discrete breathers in a Klein-Gordon chain, 
i.e. using the dynamical equations without any approxi- 
mation, a problem which has not been yet undertaken. 

In Klein-Gordon lattices, the local modes due to an 
impurity and the localized modes due to the nonlinear- 
ity of the lattice can be considered equivalent entities 
[0 |T^. This fact suggests that a moving breather can 
exhibit an interesting behaviour when interacts with an 
impurity. As a matter of fact, we have found that the 
moving breather can cross the impurity, can be reflected 
by it (being able to let the impurity excited), or can be 
trapped, generating a local depository of energy. The 
results obtained in our study show that it is necessary 
the resonance of the moving breather with a nonlinear 



local mode (i.e. a discrete static breather centred in the 
impurity site) in order to generate a trapped breather. 

This behaviour is qualitatively similar to the observed 
in the NLS approximation However, the ranges of 
the values for which the different phenomena occur are 
rather different. Furthermore, the most remarkable dif- 
ference between their approximation and our results is 
that in they obtained that the trapping is due to a reso- 
nance with a linear local mode, whereas in our case, the 
resonance must involve a nonlinear local mode. 

In addition, we have observed the non existence of 
trapping phenomena even though resonance with the 
nonlinear local mode occur. We have also observed that 
in this case the tails of the nonlinear local mode and the 
lineal local mode have different vibration patterns. As a 
consequence, we conjecture that both tails must have the 
same vibration pattern in order that the trapping occur. 



II. MODEL AND SOLUTIONS GENERATION 

A. Formulation of the model 

In order to study the effects of impurities on breather 
mobility, we consider a simple model where moving 
breathers can be generated, that is, a Klein-Gordon 
chain with nearest neighbours attractive interaction ^, 
M , whose Hamiltonian is given by: 



H 



E 



itt^ + Vn{Un) -f ^C(u„ - Un-lY 



(1) 



where u„ represents the displacement of the n-th particle 
with respect to its equilibrium position, C is a coupling 
constant and Vn{un) is the on-site potential at the n-th 
site. We choose V as the Morse potential, i.e., Vn{u) — 
Dn{e~" — 1)^, which proves to be a very suitable one to 
obtain moving breathers 0, ^ . Dn represents the well 
depth in the n-th site. An impurity will be introduced 
in the chain by means of an inhomogeneity in the 0-th 
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well, i.e., D„ = Do{l + a5„^o), with a being a parameter 
which gives account of the magnitude of the impurity. It 
takes its values in the interval [—1, oo). Hereafter, we will 
consider Do = 1/2. Notice that the impurity can be also 
introduced in the coupling. We have checked that the 
results are qualitatively the same and will be published 
elsewhere. 

This kind of chain has been used as a simple model of 
DNA dynamics; in that context it is usually referred to 
as Peyrard-Bishop model 1 13 . In the framework of this 



model, the variables u„ represents the transverse stretch- 
ing of the hydrogen bonds connecting the two bases, D 
is the dissociation energy of a base pair, and C is the 
stacking coupling constant. 

This Hamiltonian leads to the dynamical equations 

F{{Un}) = iln + + C{2Un - M„+i - U„_i) = 0. 

(2) 

These equations have two kinds of solutions, linear 
ones, which correspond to oscillations of small ampli- 
tude, and nonlinear ones, which correspond to intrinsic 
localised modes (discrete breathers). 



B. Linear modes 

The dynamical equations can be linearised if the am- 
plitude of the oscillations is small. Thus, the equations 
(H) can be written as: 

iln + ^lu„ + C(2u„ - Un+l - Un-l) = 0, (3) 

where ojn is the natural frequency of the n-th oscillator 
in the harmonic limit. It is related to the well depth of 
the Morse potential in the form = 2Z?„, which implies 
that = co'^il + aSn.o), with lOq — I sls Do has been 
chosen to be 1/2. 

Thus, this equation has — 1 solutions (being N the 
number of particles) corresponding to linear extended 
modes (LEMs) and one localised solution, which corre- 
sponds to a linear local mode (LLM) whose origin lies in 
the introduction of the impurity. 

The frequencies of the LEMs can be calculated sup- 



posing that they are plane waves (u„(t) 



„iLj{q)t-nq 



and that the LLM decays in the space following a de- 
pendence of the form Un(t) = e^'^^tj." jj^^ Thus, the 
frequencies of the LEMs are given by: 



uj{q) 



4C sin- 



2 1 



(4) 



where q is the LEM wave vector and takes its value in 
the interval (0, tt] if a < and in [0, tt) if a > 0. 

The frequency of the LLM is given by the relation jl^ : 



(5) 



It is worth remarking that LLMs with a < have wave 
vector q = TT (the particles vibrate in zigzag) while the 



FIG. 1: Frequencies of the nonlinear modes in function of 
parameter a. They correspond to C = 0.13 although the 
dependence is qualitatively similar for every value of C. 



solution for a > corresponds to q = (the particles 
vibrate in phase) . Figure |^ shows the dependence of the 
frequencies of the linear modes with a, where the isolated 
frequencies corresponds to the LLMs. 

The knowledge of the frequency of the LLM is very 
important in order to explain the properties of moving 
breathers, because discrete breathers and local modes 
can be viewed as equivalent entities when the interplay 
between nonlinearity and inhomogeneity is considered 
00- 



C. Static and moving breathers 

A static breather can be obtained solving the full dy- 
namical equations. It can be done using common meth- 
ods based on the anticontinuous limit |l^. The imple- 
mentation of these methods basically consists in calculat- 
ing the orbit of an isolated oscillator at a fixed frequency 
ujh, and using this solution as a seed for solving the 
complete dynamical equations by means of a Newton- 
Raphson continuation method. 

If the oscillator initially chosen is the corresponding to 
the impurity, a static breather centred in the impurity is 
obtained. It will be called nonlinear local mode (NLLM), 
in analogy to the linear local modes (LLMs) that appear 
in the spectrum of linear modes (equation 

Once a static breather is calculated, it can be moved 
under certain conditions. There exists a systematic 
method to calculate moving solutions ||, ||] which con- 
sists in adding to the velocities of the static breather a 
perturbation of magnitude A colinear to the direction of 
the pinning mode, and letting the system evolve in time. 
This perturbation breaks the shift translational symme- 
try of the system. In addition, the coupling of the static 
breather must be strong enough [Q. 

In this paper, we have studied breathers with a fre- 
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FIG. 2: Different regimes in the interaction of the moving 
breather and the impurity (see text). 



quency ojh = 0.8 and a coupling C = 0.13 although this 
value has been increased in order to make some compar- 
isons. These values of the parameters provide moving 
breathers with low phonon radiation for values of the 
perturbation A < 0.2. We have considered different im- 
purities given by different values of parameter a in the 
range — 1 < a < 1. 



III. INTERACTION OF MOVING BREATHERS 
WITH IMPURITIES 

A. Numerical observations 




We have studied the behaviour of moving breathers 
when they interact with an impurity. This study is per- 
formed varying the value of the magnitude of the im- 
purity a, and shows an interesting behaviour. We have 
found some critical values of parameter a (ai < a2 < 
as < 0) which separate the different regimes (see figure 

i): 

1. The impurity acts as a potential barrier. It occurs 
whenever a > and a € (—1, ai) with ai < 0. The 
breather rebounds when reaches the impurity and 
let it excited during a brief time lapse. The ampli- 
tude of this excitement decreases when the impu- 
rity is higher, i.e., when \a\ increases. Note that if 
a ~ 0, the potential barrier can be observed, i.e., 
the breather can cross the impurity provided the 
translational velocity is high enough, as it occurred 
in §. 

2. The impurity is excited and the breather rebounds. 
It occurs for a G (ai , 0:2). The energy of the excited 
impurity is higher than the energy of the NLLM 
for the considered value of a and whose frequency 
is the same as the incident moving breather, ujh. It 
implies that the excited impurity will vibrate with 
a frequency smaller than because the on-site 
potential is soft. This behavior is shown in figure 

I 

3. The breather is trapped by the impurity. It occurs 
in the interval a € (02,03). When the breather 
is near the impurity, it is absorbed and remains 



FIG. 3: Interaction of the breather with the impurity for 
a = —0.52 and A = 0.1. Top: Evolution of the moving 
breather. Bottom: Evolution of the averaged energy density 
of the impurity. En is the energy of the NLLM. 

trapped while its center oscillates between the 
neighbouring sites. The trapping of the breather 
is shown in figure |4[ Furthermore, the trapped 
breather is quasi-periodic, as can be deduced from 
its Fourier spectrum (Figure ^), and, as a conse- 
quence, it emits a great amount of phonon radia- 
tion. 

4. The impurity acts as a potential well. It occurs 
when as < a < and manifests as an accelera- 
tion of the breather when reaches a site near the 
impurity, and a deceleration when the impurity is 
crossed. 

All of these cases have been studied using a breather 
with cjb = 0.8 and C = 0.13. The translational velocity 
of the moving breathers has been chosen small so as not 
to have high phonon radiation. The critical values of 
a cannot be determined because the transition between 
the different behaviours are diffuse. An estimation of the 
critical values for = 0.8 and C — 0.13 are: ai w 
—0.54, a2 ~ —0.49 and 03 « —0.02. These regimes 
have also been found for different values of the frequency, 
although the value of 02 decreases with u^. 

If the coupling is higher, the phonon radiation is not 
so low and can mask some of the effects. 
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FIG. 4: Evolution of the moving breather for a — —0.3 and 
A = 0.1. It can be observed that the breather is absorbed by 
the impurity and gets trapped; afterwards, the breather emits 
phonon radiation and its energy centre oscillates between the 
sites adjacent to the impurity, 
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FIG. 5: Fourier spectra for the trapped breather of figure |^ 
(top) and a static breather with the same values of C and toh- 
It can be observed the existence of only one peak in the last 
case where there exist a great number of peaks in the first 
case, which is an evidence of quasi-periodicity. 



B. Justification of some results 

Some of the results exposed in the last subsection can 
be explained from the properties of the NLLMs. Con- 
cretely, if a continuation of the static breathers is per- 
formed varying the parameter a, a bifurcation appears 
for a = Oc > and another one for a = ares < 0. 
The first one is originated by a localized Floquet eigen- 
mode which abandons the unit circle. In the second 
case, the breather bifurcates with the stationary solution 
through a pitchfork (figure In this bifurcation, the 
outer branches correspond to NLLMs either with u„ > 
for t = and with u„ < for t = 0, while the central 
branch corresponds to the stationary solution, i.e. all the 
oscillators at rest. This kind of bifurcation is different 
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FIG. 6: Pitchfork bifurcation for tOb = 0.8 and C = 0.13. 
ares = -0.5328. 



to the broken pitchforks obtained in disordered systems 
[ prl |lj, |l^. It is due to the fact that, even thought an 
isolated impurity breaks the shift translational symme- 
try of the system, it does not break the mirror symmetry. 
Consequently, the pitchforks are not broken. The value 
ares coincides with the point at which the frequency of 
a LLM is the same as the frequency of the NLLM, i.e., 
wl = Wb- This value can be calculated from the equation 
^ as a function of Wb and C : 

«re. = -^H^u;l){u;l~u;l~4C) (6) 

Thus, for C = 0.13 and uj^ = 0.8, ares = -0.5628 
which is lower than ai. We have also checked for different 
values of the frequency and coupling and we have found 
that there are situations where ai is slightly smaller 
that ares but always Ures < a2- However, the trapped 
breather does not exist for a > 0. It indicates that 
a G {ares,0) is a necessary condition for the existence 
of the trapped breather. 

A conjecture for the existence of trapped breathers 
when a < is the following: in this case, the LLM 
has g = 0, and also all the particles of the NLLM vi- 
brates in phase; this vibration pattern indicates that the 
NLLM bifurcates from plane waves with q = |l^, i.e., 
the NLLM bifurcates from the LLM and it will be the 
only localized mode that exists when the impurity is 
excited for a > ares- In fact, we have performed suc- 
cessful a continuation from the NLLM to the LLM at 
constant action and disorder |l^, |l^. Thus, when the 
moving breather reaches the impurity, it can excite the 
NLLM. However, Doppler effect makes the NLLM 'see' 
the moving breather with a frequency different from uj^ 
so that when they interact, they merge into an entity 
(the trapped breather) with two frequencies. 

When ares < a < a2, the NLLM is unable to create 
a trapped entity. It is observed that the energy of the 
NLLM decreases with so, there must be threshold 
for the energy of the NLLM in order that the trapped 
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breather can exist. The narrow window of impurity ex- 
citements observed in the interval (01,02) can be due 
to a resonance of the moving breather with the NLLM 
of frequency shghtly smaller that and, therefore, the 
excited impurity has a higher energy than the NLLM of 
frequency ujh, as we are considering a soft potential. 

For a < ares, the trapped breather cannot be gener- 
ated, and the moving breather always rebounds. In addi- 
tion, the NLLM does not exist. Therefore, there must be 
connection between both facts, i.e., the existence of the 
NLLM will be a necessary condition in order to obtain a 
trapped breather. 

If a > 0, the scenario is different. In this case, the 
LLM has q = tt but the NLLM's particle vibrate again 
in phase, that is, the NLLM will not bifurcate from the 
LLM. Thus, there are two different localized modes for 
a > 0: the LLM and the NLLM. But, actually, the equa- 
tions that govern the system are nonlinear, so the linear 
modes can only correspond to low-amplitude oscillations. 
In the case of the NLLM, the linear regime corresponds 
to the tails. Thus, if the moving breather reaches the 
impurity site, it will excite the NLLM, its tails and also 
the tails of the LLM. But the latter vibrate in zigzag. As 
a consequence, there will be two different linear localized 
entities: the tails of the LLM (vibrating in zigzag) and 
the tails of the NLLM (vibrating in phase). Therefore, 
we conjecture that the existence of both linear localized 
entities at the same time may be the reason why the 
impurity is unable to trap the breather when a > 0. 

To summarize, the existence of a NLLM for a cer- 
tain value of a is necessary for the existence of trapped 
breathers. However, if there exists a LLM with a differ- 
ent vibration pattern to the NLLM, the trapped breather 
does not exist. 



IV. CONCLUSIONS 

We have studied the interaction of a moving discrete 
breather in a Klein-Gordon chain with an impurity at 
rest. A rich behaviour is observed when the breather 
reaches the impurity. It can be summarized as follows: 
a) the impurity can act as a potential barrier or a po- 
tential well; b) the breather can be reflected by the im- 
purity which is let excited; c) the impurity can trap the 
breather, making it quasi-periodic. 

A previous work made using the NLS approximation 



pOl showed that moving breathers interacting with an 
impurity have qualitatively the same behaviour as the 
observed by us. There is, however, a very important 
difference between this work and ours. It relies in the 
fact that in the NLS approximation, a resonance between 
the frequencies of the moving breather and the linear 
local mode created by the impurity is needed whereas 
in our case, this resonance prevents the breather from 
being trapped. In our case, the trapping phenomena is 
owed to the resonance between the moving breather and 
a nonlinear localised mode. 

The model used in our study is the same as the stud- 
ied by Peyrard and Bishop to explain DNA denaturation 
[ p^ . It indicates that our results can be applied to study 
some properties of DNA chains. For instance, if a moving 
breather is generated in a DNA chain, it can be trapped 
by an impurity, and act as a precursor of the transcrip- 
tion bubble This trapped breathers can interact 
with another breathers and collect energy |p^ . 

The impurity in DNA can have two different origins. 
The first one consists in a modification of the well depth 
due to the action of a transcription enzyme through a 
chemical effect, as explained in |T^. The second one 
relies in the fact that the A-T base pairs have 2 hydro- 
gen bonds whereas the C-G base pairs have 3 hydrogen 
bonds. So, it can be supposed that in the first case, the 
well depth of the Morse potential is 2/3 of the second 
case. From the results of this paper, trapping does not 
occur when the well depth in the impurity is higher than 
in the homogeneous area. It implies that trapping can 
occur in a chain of C-G pairs with an impurity of A-T. 
Nevertheless, the main fault of these mechanisms is the 
assumption of homogeneity in DNA, where the role of 
the inhomogeneity due to the genetic code is known to 
be crucial in the dynamics ||l9|. However, the study of 
the effects of impurities is the first step to understand 
the dynamics of moving breathers in DNA sequences. 
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